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ABSTRACT 

We discuss two multiple star systems that host known exoplanets: HD 2638 
and 30 Ari B. Adaptive optics imagery revealed an additional stellar companion 
to both stars. We collected multi-epoch images of the systems with Robo-AO 
and the PALM-3000 adaptive optics systems at Palomar Observatory and provide 
relative photometry and astrometry. The astrometry indicates that the compan¬ 
ions share common proper motion with their respective primaries. Both of the 
new companions have projected separations less than 30 AU from the exoplanet 
host star. Using the projected separations to compute orbital periods of the new 
stellar companions, HD 2638 has a period of 130 yrs and 30 Ari B has a period of 
80 years. Previous studies have shown that the true period is most likely within 
a factor of three of these estimated values. The additional component to the 30 
Ari makes it the second confirmed quadruple system known to host an exoplanet. 
HD 2638 hosts a hot Jupiter and the discovery of a new companion strengthens 
the connection between hot Jupiters and binary stars. We place the systems on 
a color-magnitude diagram and derive masses for the companions which turn out 
to be roughly 0.5 solar mass stars. 

Subject headings: binaries: visual - instrumentation: adaptive optics - stars: individual 
(HD 2638, 30 Ari B) stars: solar-type 
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1. INTRODUCTION 

Since the discovery of the first exoplanet 
system in the 1990s, planetary systems 
have been found in an increasingly varied 
array of architectures. With this variety 
has come an increasing number of mod¬ 
els and theories to explain these systems. 
It has also been realized that stellar com¬ 
panions play a key role in the evolution of 
the planetary dynamics for some systems. 
There are several proposed ways in which 
binary stars can influence the orbital prop¬ 
erties of exoplanets. 


One of these is Kozai migration ( 

Wu & Murrav 

2003; 

Wu et al. 2007: 

Fabrvckv & Tremaine 


2007; 

Takeda et al. 

2008) where the mu- 



tual torques between the binary and the 
planet transfer angular momentum, caus¬ 
ing the orbit of the planet to become more 
eccentric. Eventually the eccentricity be¬ 
comes high enough that the planet ap¬ 
proaches the star, which then tidally cir¬ 
cularizes the orbit. Not all exoplanets in 
binary systems undergo Kozai migration, 
as it depends on the mutual inclination 
of the orbits as well as the architec t ure o f 
the planetary system. Naoz et al. ( 2012 ) 


found that it can account for about 30% 
of the observed hot Jupiter planets, which 
matches well with the projected spin-orbit 
angle distribution of hot Jupiters. 

The stellar companion can also inter- 
act w ith the protoplanetary disk. Rahkov 


(2Q13) proposed that the stellar compan¬ 
ion effects the protoplanetary disk by slow¬ 
ing the planetesimals in the disk and allow¬ 
ing the formation of planets. This requires 
a relatively massive disk, but can produce 
syste ms such as 7 Cep and a Cen ([Rahkov 

201 oh . 


Kiev fe Nelson] (2008) suggest an alter¬ 
native method in which an eccentric stel¬ 
lar companion leads to a strong periodic 
disturbance in the disk whenever the com¬ 
panion is at periastron. This leads to an 
eccentric disk that changes over time due 
to the interactions with the stellar com¬ 
panion. This can produce massive eccen¬ 
tric planets, depending on the initial posi¬ 
tion of the planetesimal. Their modeling 
did not consider the initial creation of the 
planetesimals. 

Circumbinary planets have also been 
discovered, where the planet orbits both 
stars in a shor t-period orbit, on the order 
of tens of days ( Lee et ahll2009l ; Dovlc ct al. 


2 0111 ). In these systems the planet (s) and 
the stellar companion are all co-planar 
implying that they formed from flat cir¬ 
cumbinary disks. Since their discovery 
there is considerable effort to model the 
formation and dynamics of th ese syste ms 
( Kiev fe Haghighipourll2014 ; Armstrong et al. 


20141 ) . 


In the case of multiple stellar systems, 
coplanar systems are rare. Of the 2413 sys¬ 
tems cataloged in the Sixth Catalog of Or- 


bitsof Visual Binary Stars (lHartkopf et al. 


2001 ) 48 systems have orbit solutions with 


all seven Campbell elements for both inner 
and outer systems in hierarchical arrange¬ 
ments. Of these, only four have the pos¬ 
sibility of being coplanar according to the 
precepts of lFekell ( 198lh . 

The various models of exoplanet forma¬ 
tion in binary systems can only be fully 
tested when there are sufficient numbers 
of known systems, especially systems with 
relatively short binary periods. There have 
been a number of surveys to measure the 
duplicity rates among the exoplanet hosts 
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(Patience et al. 

2002 : 

Eggenberger et al. 

2007 


Ohauvin et al 

20 ll[ 

Mason et al. 

2011 


Roberts et al. 

2011 ; 

Ginski et al. 


2012). These surveys have revealed sev¬ 


eral interesting features. 

The frequencies of exoplanets among 
single stars and components of wide bi¬ 
naries (semi-major axis larger than 100 


AU) are indistinguishable (I 

laghavan et al. 

2006: Bonavita & Desidera 

2007 

). Also, 


the population of exoplanets in wide bi¬ 
naries is essentially the same as the popu¬ 
lation around single stars, indicating that 
wider binaries have minimal impact on dy¬ 


namics of the exoplanets OBonavita fc Desidera 


2007). The population of exoplanets in bi¬ 


naries with semi-major axis smaller than 


orbiting single stars (7 

bicker & Mazch 

2002 ; 

Bonavita & Desidera 

2007). There 

are 


fewer exoplanets in tight stellar binary 
systems, but the exop la nets te nd to be 


more massive (iDesidera fe Barbieril 120071: 


Duchen e 201 y). No planet has been de¬ 


tected in stellar binaries that have a sep¬ 


arati on of less than 10 AU flRoell et al. 


2012). That same paper also determined 
that multi-exoplanet systems have only 
been detected in binary stars with a pro¬ 
jected separation larger than 100AU. 

Known close visual binaries are tra¬ 
ditionally excluded from radial-velocity 
(RV) exoplanet programs because the pres¬ 
ence of a visual companion degrades the 
RV precision. This intrinsic bias compli¬ 
cates statistical inferences about exoplan¬ 
ets in binaries. Moreover, a faint visual 
companion that is itself a close (spectro¬ 
scopic) binary pair can produce periodic 
low-amplitude RV modulation in the com¬ 
bined light that can be mistaken for an 


exoplanet. False positive exoplanet detec¬ 
tions caused by unrecognized hierarchical 
multiplicit y of their hosts m ay reach 1— 
2 percent ( Tokovinfnl 2014b 1 of the total 
exoplanet sample. This is yet another rea¬ 
son to observe exo-hosts with high angular 
resolution and deep dynamical range. 

In addition to dynamical interactions, 
exoplanets and their host stars share a 
common origin, and determining the prop¬ 
erties of that local environment for the star 
leads to information about the exoplanet. 
Examples of this include metallicit y co r 
relations (e.g. 16 Cyg AB, ISchuler et al. 


201 l|), stellar rotation, and activity. 

This paper details two multiple star sys¬ 
tems that host known exoplanets: HD 2638 
and 30 Ari B. Additional stellar compo¬ 
ne nts were found in ea ch of these systems 
by Rid dle et_ ah (j2015l) . This paper dis¬ 
cusses the systems and those newly discov¬ 
ered components systems in depth. The 
observations of the systems are covered in 
Sectional Section [3] discusses the HD 2638 
system and Section [4] discusses the 30 Ari 
system. Our results are summarized in 
Section 0 

2. OBSERVATIONS 
2.1. ROBO-AO 

The companions to HD 2638 and 30 Ari 
B were first detected using the Robo-AO 

system as part of a survey of nearby bina : _ 

ries w ith solar type components (IRiddle et al 
20151 1. We have copied the astrometry and 
photometry from that paper in order to 
do a more in-depth analysis of these two 
stars. For those unfamiliar with Robo- 
AO, we provide a brief introduction to 
the system. Robo-AO is a robotic laser 
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guide star adaptive optics (AO) system 
( Baranec et al.1 2013 . 2014 ) that operates 
automatically in order to carry out high ef¬ 
ficiency observing programs. Robo-AO is 
currently deployed on the Palomar Obser¬ 
vatory 60-inch telescope. Robo-AO uses 
Rayleigh scattering from a UV laser fo¬ 
cused at 10 km from the telescope as the 
wavefront reference, and generates diffrac¬ 
tion limited images with Strehl ratios of 4- 
26% in the i filter. The AO system corrects 
the high order wavefront aberrations with 
automated software that operates at a rate 
of 1.2 kHz, sharpening the instantaneous 
point spread function across the science 
camera field of view (FOV). A bright star 
within the FOV is still required to correct 
the tip-tilt motion, which is corrected by 
the automated data processing software. 

The Robo-AO science camera is an elec¬ 
tron multiplying CCD (EMCCD) detec¬ 
tor with 1024x1024 pixels; the pixel scale 
is 43.1 mas/pixel, with a total FOV of 
44"x44". The EMCCD is operated at 8.6 
frames per second in science mode, with 
an exposure time of 0.115 s each. A total 
of 516 frames are gathered during each 60 
s exposure, which are then combined into 
a single image for further analysis by the 
automated data processing software. Both 
tar gets were observe d with the SDSS i 'fil¬ 
ter ( York et al.12000 ). and additionally HD 
2638 was observed with the SDSS r' and 
z' bands. The data were collected dur¬ 
ing multiple observing runs in 2012 and 
2013. See Sections [3] and Q] for details. 
The Robo-AO data are procesed by the 
stan dard pipeline, which is discussed in 


if 

Riddle et al.l fl2015l ). That paper did not 


publish error bars for the astrometry mea¬ 
sured using speckle processing. It did give 


mean errors on measurements of control bi¬ 
naries. These companions are fainter than 
the control binaries, so to be conservative 
we have set the error bars as 10 mas, which 
translates to 1?1 in position angle at 0'.'5 
separation. 


2.2. PALM-3000 Adaptive Optics 


In order to confirm and further charac¬ 
terize the companions, we observed both 
stars on 2013 September 28 UT with 
the Palomar Observatory Hale 5 m tele¬ 
scope using the PALM-3000 AO system 
and the PHARO near-IR camera. The 
PALM-3000 AO system is a natural guide 
star system using tw o deformable mirrors 
f Dekanv et al.1 2013 ;). One corrects low- 
amplitude high spatial frequency aberra¬ 
tions, while the other corrects the higher- 
amplitude low spatial frequency aberra¬ 
tions. It is optimized for high contrast ob¬ 
servations and routinely produces Strehl 
ratios greater than 80% in Ks band. The 
PHARO camera uses a HgCdTe HAWAII 
detector for observations between 1 and 


2.5/mi wavelength (J Iavward et al. 2001). 
The camera has multiple filters in two fil¬ 
ter wheels. For HD 2683 we collected 100 
frames in both the J and Ks filters, while 
we collected 50 frames of 30 Ari B in the 
same filters. After debiasing, flat fielding, 
bad pixel correction and background sub¬ 
traction, the frames were coadded. After 
the individual frames were calibrated, we 
created multiple images, by coadding 10 
frames into each image. The iterative blind 
deconvolution algorithm, fitstars, was used 
to measure the astrom etry and photome¬ 


try of the objects flten Brummelaar et al. 
1006. 2 ()()()!) ■ 


Photometric error bars were assigned 
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using the technique described in lRoberts et al. 
( 20051 ). For the astrometric error analy¬ 
sis, we computed the standard deviation 
of the astrometric measurements. Analy¬ 
sis of simula ted binaries created u sing the 
technique of iRoberts et al.l (120051 ) showed 
that this error is almost entirely measure¬ 
ment error. The other major sources of 
astrometric error are the errors in the cali¬ 
bration of the plate scale and the position 
angle offset. The PHARO pixel scale has 
been measured many times with a combi¬ 
nation of an astrom etric mask and binary 
star measurements flMetchevI l2006h . The 
changes between observing runs are small. 

On the night these data were acquired, 
there was not enough time to repeat the 
lengthly astrometric mask calibration pro¬ 
cess, but we observed six calibration bina¬ 
ries and updated the plate scale and the 
orientation angle offset. This resulted in 
a plate scale of 0"025±0"0005 and a po¬ 
sition angle offset of 0?7A0?5. The final 
error bar on the astrometry was the root 
sum square (RSS) of the calibration errors 
and the scatter in the astrometric mea¬ 
surements. 

3. HD 2638 


HD 2638 (HIP 2350 = WDS 00293- 
0555) is a nearby sola r type star with an es - 
timated age of 3 Gyr (IMoutou et al.ll20051) . 

It is paired with HD 2567 (HIP 2292), 
another solar type star, in a wi de (8 39") 
comm on-proper-motion binary. IShava fe Oiling 
(20111) determined that the pair had a 99% 
chance of being true companions through a 
Bayesian analysis of the measured astrom¬ 
etry. The astronomical details of the HD 
2638 system are given in Table [T] The mea¬ 
sured parallaxes from Hipparcos of the two 


stars are slightly different. For HD 2638 
it is 20.03±1.49 mas and for HD 2567 it 


is 17.63i0.79 mas (Ivan Leeuwenl 120071) . If 
the parallaxes are correct, it would mean 
that the stars are separated by 6.8 pc, 
which would mean they are not physi¬ 
cally bound. This woulcL disagree with 
the conclusion of Shava & Oiling! fcoill ). 
Independent confirmation that A and B 
are physical comes from their measured 
RVs being almost identical. A9.60A0.3 


km/s for A ([Nordstrom et al.l 12004), an d 


A9.61i0.01 for B (jLatham et al. 


2002 ). 


For the distance to the system, we will 
use the value for HD 2567 since it has the 
lower error bar. Unresolved binary stars, 
such as turns out to be the case with HD 
2638 (see below), can cause errors in par¬ 
allaxes and proper motions measured by 
Hipparcos due to o rbital motion in th e case 
of close pairs ( Shatskii fc Tokovinin 19981) 
and due to a cross- talk with the scannin g 
law for wider pairs (Tokovinin et al. 20131) . 
This most likely explains the difference. 

Radial velocity (RV) observations de¬ 
tected a hot Jupite r orbiting HD 2638 


(IMoutou et al.l 120051) . The planet, HD 


2638b, has a minimum mass of 0.48Mj up 
and an orbital period of 3.4442 i 0.0002 
d. The orbital semi-major axis is 0.044 
AU and the orbit has been tidally circu¬ 
larized. The planet’s equilibrium temper¬ 
ature is on the order of 1100 K (with an 
assumed albedo of 0.3). 


Riddle et ah (120151 ) imaged a compan¬ 


ion to HD 2638 on 2012 September 3 UT. 
Over the next 13 months, two more ob¬ 
servations were carried out with Robo-AO 


Riddle et al.l 020151 ) and another observa¬ 
tion was made with the 5 m Hale tele¬ 
scope using the Palm 3000 AO system and 
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Table 1 

Basic information on two multiple systems 


Object 

WDS 

(J2000) 

HD 

HIP 

TfflP 

(mas) 

V 

(mag) 

Spectral 

type 

HD 2567 A 

00293-0555 

2567 

2292 

17.63+0.79 

7.76 

GO 

HD 2638 B 


2638 

2350 

20.02+1.49 

9.37 

G5 

30 Ari A 

02370+2439 

16246 

12189 

23.93+0.59 

6.48 

F5V 

30 Ari B 


16232 

12184 

24.52+0.68 

7.09 

F6V 


the near-IR PHARO camera. Images of 
HD 2638 BC are presented in Figure [21 
The measured astrometry and photome¬ 


try ar e shown in Table [2j In iRiddle et al. 


020151) there are two measurements from 
2012.7621. One has a position angle of 
167?3 and the other is 170?0. The later 
measurement clearly is in error, as it does 
not fit the progression of the astrometry. 
We have chosen to ignore that astrometric 
point. 

HD 2638 has a proper motion o f 0.247 
± 0.0012"/yr ( van Leeuwen 2007 1. We 
compare this to the relative motion of the 
companion to the primary of 0.0108 ± 
0.0175"/yr. It is also in a different di¬ 
rection as shown in Figure CD While the 
error bar on the relative motion is high 
compared to the motion, both are much 
smaller than the proper motion and we 
conclude that the companion has common 
proper motion with HD 2638. The mea¬ 
sured motion is consistent with orbital mo¬ 
tion in that both the position angle and 
separation are changing in a regular man¬ 
ner within the error bars. With only four 
data points points it is impossible to de¬ 
termine if there is a deviation from linear 
motion, which would be required to con¬ 


firm orbital motion. The binary pair is 
designated HD2638 BC. Normally the pair 
would be designated HD 2638 Ba,Bb, but 
that could cause confusion with the plane¬ 
tary companion. 


3.1. DISCUSSION 


The components of the newly discovered 
binary system can be placed on the color- 
magnitude diagram (CMD), using the pho¬ 
tometry in the combined light from the lit¬ 
erature and the differential photometry we 
measured with Robo-AO and PHARO. We 
also use the the Dartmouth isochrones for 
1 G yr and 3 Gyr age w ith solar metallic- 
ity ( Potter et ah 2008)0. The isochrones 


define the relations between various col¬ 
ors and between mass and absolute mag¬ 
nitude. The combined J and K magni¬ 
tudes are taken from 2MASS, but the com¬ 
bined SDSS i magnitudes for our targets 
are not measured directly. We determined 
i by extrapolating from the known V and J 
magnitudes, using polynomial relations be¬ 
tween the V — i, i — J, and V — K colors de¬ 
rived from the Dartmouth isochrones. The 
combined, differential, and individual mag- 


1 http://stellar.dartmouth.edu/models/webtools.html 
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Fig. 1.— The measured astrometry for HD 2638 BC is plotted. The black arrow is the 
position the companion would have been at our last observation if it was a background star 
with no proper motion. The red line shows the actual motion between our first and last 
observation. 


nitudes of the components are listed in Ta¬ 
ble [3] 


Figure [3] shows the position of the bi¬ 
nary components in the (z, i — J) CMD, 
together with the 3-Gyr isochrone defining 
the main sequence. The masses in Table [3] 
are determined from the isochrone and the 
absolute magnitudes in the z band. 


We estimate the orbital period, P, of 
the binary from the measured separations, 
p, distances, and mass sum A4 using Ke- 

3 

pier’s Third Law, = M.. The median 
ratio between the projected separation, d 
(= —2—), and orbital semi-maior axis, a, is 
close to one flTokovinin 2014af ). with scat¬ 
ter by a factor of two caused by orbital 
phase, orbit orientation, and eccentricity. 
The strict lower limit is a > |. Statisti¬ 
cal period estimates using the assumption 
a = d are denoted as P*. 


The masses of HD 2638 B and C esti¬ 
mated from their absolute i! magnitudes 
correspond to the spectral types G8V and 


M1V. The projected separation of 28.5 AU 
leads to P* ~ 130 yr. The true period 
likely equals P* to within a factor of three. 
There is however a small probability that 
the true separation is much larger than 
the projected separation, in which case the 
actual period can be substantially longer 
than P*. 


The projected separation of HD 2638 
B and C is quite close. The stellar com¬ 
panion has most likely had a significant 
impact on the dynamical evolution of the 
planet. It has long been thought that 
hot Jupiters formed beyond the snow line 
and migrated to their current locations. 
There are a number of theories on what 
caused that migration. One suggested 
mechanism is that interaction between the 
protoplanetary disk and the plan et causes 


the planet to migrate inwards (Lin ct ah 


19961 ). A second method is planet-planet 


scattering that puts exoplanets into ec¬ 
centric orbits that are then tidally cir- 
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Table 2 

Relative astrometry and photometry of HD 2638 BC. 


Epoch 

on 

pH 

A r 

A i 

Az 

A J 

AKs 

Instrument 

2012.6748 

166.4T1.1 

0.528T0.010 


3.38 




Robo-AO 

2012.7621 

168.7T1.1 

0.553T0.010 

3.87 


3.06 



Robo-AO 

2013.6208 

166.7T1.1 

0.526T0.019 


3.19 




Robo-AO 

2013.7422 

167.3T1.4 

0.52 ±0.01 




2.50T0.09 

1.99T0.05 

PHARO 


Table 3 

Photometry of the new pairs 



i 

(mag) 

J 

(mag) 

Ks 

(mag) 

Mass 

(M e ) 

HD 2638 BC, m 

— M = 

= 3.77 mag 


BC combined 

8.91 

6.08 

5.82 


BC Am 

3.29 

2.50 

1.99 


B 

8.96 

7.96 

7.47 

0.87 

C 

12.25 

10.46 

9.46 

0.46 

30 Ari 

BC, rn - 

■ M = 

3.05 mag 


BC combined 

6.88 

6.08 

5.82 


BC Am 

4.31 

3.44 

2.84 


B 

6.90 

6.13 

5.90 

1.11 

C 

11.20 

9.57 

8.74 

0.50 
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Fig. 2.— On the left is an i image of the HD 2638 BC from Robo-AO. It is a subimage of 
the full 44" held of view and is ~2"5 across. An arrow points to the location of the newly 
discovered companion. On the right is a Ks image of the HD 2638 BC acquired with the 
PALM 3000 AO system. It is a subimage of the full 25" held of view and is ~ 2" across. 


cularized flC hatte rjee et al. 2008]). An¬ 
other method is Kozai migration, where 
a stellar companion on an inclined orbit 
induces exoplanets into eccentric orbits. 
The primary star will then tidally circu ¬ 
larize the planet f Wu fe Murray 2003 ). 
Nagasawa et al.l ( 2008 ) suggested that 
planet-planet scattering could lead to ec¬ 
centric inclined exoplanets, which would 
then start Kozai cycles with the other ex¬ 
oplanets in the syst em resulting in hot 
Jupiters. Wu fe Lithwick ( 201 ill proposed 
that secular chaos could cause the inner¬ 
most planet to lose angular momentum 
and then be tidally circularized forming a 
hot Jupiter. It is quite probable that all 
of these methods occur in nature, but we 
do not know which is the dominant mech¬ 


anism. 


The presence of a stellar companion to 
HD 2638 does not help us differentiate be¬ 
tween these methods. While Kozai migra¬ 
tion requires a third body in the system, 


the other methods are not ruled out by the 
existence of a stellar companion. To gain 
a deeper understanding of the dynamical 
evolution of this system we will need to 
monitor the system to eventually compute 
an orbit of the stellar system. Also con¬ 
tinued RV monitoring is needed to identify 
any additional planets in long periods or 
highly inclined orbits. Continued monitor¬ 
ing would allow the measurement of the 
RV signature of the companion star which 
will allow a more complete orbit to be de¬ 
termined. 

It is conceivable that the newly discov¬ 
ered component, C, is a close binary mas¬ 
querading as an exo-planet. We estimate 
that C is fainter than B by about 4 mag 
in the visible, so the observed RV ampli¬ 
tude of 67.4 m s -1 in the combined light 
of B+C corresponds to the 2.7km s^ 1 am¬ 
plitude in the velocity of C. The minimum 
mass of the hypothetic binary companion 
Cb is then 0.012 solar, so it could be a 0.1 
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Fig. 3.— Position of the binary components in the color-magnitude diagram. The full line 
shows a 1 Gyr isochrone (correspondong to 30 Ari B) and the dashe_d line : shows the 3 Gyr 
isochrone (corresponding to HD 2638), both with solar metallicity flDotter et al.l 20 081 ). 


A'f© dwarf in a 3.44-day orbit with a 30° 
inclination. This is possible, although sta¬ 
tistically unlikely (probability 0.13). If, on 
the other hand, Ca and Cb had comparable 
masses, the RV amplitude of this putative 
binary would be about 70 km s -1 , and an 
orbital inclination of ~ 20° would suffice 
for blending with the lines of B and pro¬ 
ducing the fake exo-planet RV signal. In 
the light of the discovery of C, a more in- 
depth analysis of the spectra of HD 2638 
is desirable to validate the exo-planet. 

4. 30 ARI B 

30 Ari B (HIP 12184 = HD 16232 
= WDS 02370+2439) is a star classified 


as F6V f Guenther et al. 20091 at a dis¬ 


tance of 40.8+ pc (Ivan Leeuwenl 120071) . 
Along with HIP 12189, it forms the bi¬ 
nary system 30 Ari. The pair has a 
separation of ~40" and an estimate d 


perio d of 34,000 year s (ITokovinin et al. 


2006). iShaya fe Ollingl (20111 ) determined 


that the pair had a near 100% probabil¬ 
ity of being a physical system. 30 Ari 
A (HIP 12189 = HD 16246) has been 
lis ted as an F6III, but Guenther et al. 


( 20091 ) analyzed its colors and concluded 
it is actually an F5V star. 30 Ari A is 
also a spectroscopic binary with a period 
of 1.109526 d in a nearly circular orbit 
( Morbev fe Brosterhusl 1974 b 30 Ari A 

has a measured age of 0.86+0.63 Gyr and 
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30 Ari B has a measured age of 0.91 ±0.83 
Gyr (Guenther et al. 2009). Both st ars a re 


slightly metal rich 1 Guenther et al. 2009 ). 
The details of the system are also given in 
Table [H 


Using the RV technique, Guenther et al. 


( 120091 ) discovered a 9.88±0.94 M Jup exo¬ 
planet orbiting 30 Ari B with a period 
of 335.1±2.5 d and an eccentricity of 
0.289±0.092. Based on the star’s mea¬ 
sured rotational velocity and the aver¬ 
age r otati onal velocity of F-type stars, 
Guenther et all (120091) concluded that the 
star is being observed nearly equator on. 
Assuming that the planet’s orbital axis 
aligns with the star’s spin axis, the mini¬ 
mum mass is close to the true mass. With 
just over six years of data, there are no 
indications of additional signals in the RV 
data coming from additional companions. 


Riddle et al. ( 2015! ) detected a compan¬ 


ion to 30 Ari B with Robo-AO on 2012 
September 3 UT at a separation of just 
over half an arcsecond. The compainion 
was confirmed with Robo-AO approxi¬ 
mately a year later and we conducted fol¬ 
low up observations in the near-IR with 
the Palm 3000 AO system on the Palomar 
5 m telescope at approximately the same 
time. Images of 30 Ari BC are shown in 
Figure [51 The measured astrometry and 
photometry are in Table HJ The error bars 
of 30 Ari BC are smaller than the error 
bars of HD 2638 BC because 30 Ari BC 
is a brighter star and AO performance is 
improved on brighter targets. 

30 Ari B ha s a proper motion o f 0.151 
±0.00075"/yr ( van Leeuwen 2007 ). Be¬ 
tween the first and last observation the 
companion moved 0.0229 ± 0.0199"/yr. 
From this we conclude that the compan¬ 


ion has common proper motion with 30 Ari 
B. Like HD 2638, this would normally be 
designated 30 Ari Bb, but to avoid con¬ 
fusion with the planet we designate it as 
30 Ari C. The measured motion is consis¬ 
tent with orbital motion in that both the 
position angle and separation are chang¬ 
ing in a regular manner. With only three 
data points it is not possible to determine 
if there is a deviation from linear motion, 
which would be required to confirm orbital 
motion. 


As part of a duplicity survey of exo¬ 
planet host stars, Gi nski et al. (2012) ob¬ 
served 30 Ari B with the 2.2 m Calar Alto 
telescope and the AstraLux lucky imaging 
camera in the SDSS i filter on 23 February 
2010. They did not detect the companion 
due to the limited dynamic range of the 
lucky imaging system and poor weather. 


4.1. DISCUSSION 


The projected separation between B 
and C is 22.3 AU and corresponds to the 
statistical period P* = 80 yr. Using the 
same technique as used for HD 2638, we 
computed photometry and masses for 30 
Ari BC which are presented in Table El 
The luminosities and masses of B and 
C correspond to spectral types F8V and 
M1V, respectively. 

As in the case of HD 2638 B, the possi¬ 
bility that the exoplanet signal comes from 
a hierarchical multiple system n eeds to be 
addressed. According to ( Guenther et al. 


20091 ). the axial rotation of B is 38 km s -1 , 
while the RV amplitude is 0.27 km s _1 . Al¬ 
though at visible wavelengths C is fainter 
than B by at least 5 mag, it has a lower 
effective temperature and stronger lines, 
making it difficult to evaluate the ampli- 
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Fig. 4.— The measured astrometry for 30 Ari BC is plotted. The black arrow is the position 
the companion would have been at our last observation if it was a background star with no 
proper motion. The red line shows the actual motion between our first and last observation. 


tude of the hypothetic binary Ca,Cb re¬ 
quired to mimick the exo-planet. An am¬ 
plitude of 10 km s _1 , for example, can be 
produced by a companion Cb of 0.26 A4 0 
at 90° inclinantion. 


Guenther et ah f 2009h points out that, 


“...according to Eg genberger e t al. (120071 ) 
most of the very massive planets are in bi¬ 
nary systems. But ... the distance to 30 
Ari A is presumably too large to have any 
effect on the properties of the planet.” Our 
discovery of 30 Ari C may explain why it 
is not the exception. 


The detection of a stellar companion to 
30 Ari B means that the 30 Ari system 
is a quadruple system containing at least 
one exoplanet. This makes it the second 
confirmed quadruple system known to host 
an exoplanet. The first is KIC 4862625 


hosting the planet Phlb (Schwamb et al. 


2013). Phlb is a circumbinary planet or¬ 


biting an eclipsing binary consisting of G- 
and M-type dwarfs. The other two com¬ 
ponents of the quadruple are in a tight vi¬ 


sual binary with a projected separation of 
1000AU from the G/M pair. 

The similarity between KIC 4862625 
and 30 Ari is that both systems are 
quadruples consisting of two relatively 
close pairs that are widely separated. 
These 2+2 systems are relatively common. 


Riddle et al.l (1201511 found that over 10% 


of wide binaries are actually 2+2 quadru¬ 
ples and that the majority of quadruples 
are in this form, rather than the form of 
a triple system with a distant companion 
(i.e. 3+1). In both of the quadruple sys¬ 
tems hosting exoplanets, the two binaries 
making up the quadruple are most likely 
far enough apart that there is no impact 
on the dynamics or the formation of the 
exoplanet. 


5. SUMMARY 

We have detected stellar companions to 
the exoplanet host stars HD 2638 and 30 
Ari B. Both systems are part of wider mul¬ 
tiple star systems. In the case of 30 Ari B, 
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Fig. 5.— On the left is an i image of the 30 Ari BC from Robo-AO. It is a subimage of the 
full 44" held of view and is ~2"5 across. The location of the newly discovered companion 
is indicated by the arrow. On the right is a Ixs image of the 30 Ari BC acquired with the 
PALM 3000 AO system. The object to the lower right of the primary star of the image is a 
ghost image caused by the neutral density filter. This is a subimage of the full 25" held of 
view and is ~2" across. 


Table 4 

Relative astrometry and photometry of 30 Ari B. 


Epoch 

»(°) 

P(") 

A i 

A J 

A Ks 

Instrument 

2012.6750 

2013.6167 

2013.7426 

285.9 ±1.1 

285.2 ±0.2 

284.77±0.79 

0.528±0.010 

0.536±0.002 

0.55 ±0.008 

4.32±0.1 

4.31±0.1 

3.44±0.15 

2.84±0.07 

Robo-AO 

Robo-AO 

PHARO 


the other component, 30 Ari A, is also a 
binary star making it a quadruple system. 
This is the second quadruple system known 
to host an exoplanet. 


Bonavita & Desidera ( 2 007) found that 


binaries with exoplanets tend to have low 
mass stellar companions. Both of the sys¬ 
tems detected here continue that trend. 
This could be a selection effect, as most 
RV surveys have eliminated known bi¬ 
naries from their observing lists. These 


known binaries were mostly discovered by 
the Hipparcos satellite, visual microme¬ 
ters or speckle interferometry instruments, 
which all have limited dynamic range. AO 
is currently being used for most exoplanet 
duplicity searches and it can achieve a 
much higher dynamic range allowing for 
the detection of late-type stars. 

HD 2638 BC and 30 Ari BC both have 
projected separations of less than 30 AU 
and add to the rather small number of sys- 
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terns with such low separations. The stel¬ 
lar companions have the potential to have 
been a key part of the dynamical evolu¬ 
tion of the exoplanets. However, the verac¬ 
ity of these exoplanets has to be addressed 
in light of the new companions, to elimi¬ 
nate the possibility that these RV discov¬ 
eries are false positives caused by hierarchi¬ 
cal multiplicity. These systems strengthen 
the link between massive planets and bi¬ 
nary stars. The astrometry of these bi¬ 
naries should continued to be monitored. 
In the near future this will hopefully de¬ 
tect orbital motion and provide an indica¬ 
tion to the period of the objects. With a 
long enough span of observations, an or¬ 
bit can be computed, which will provide a 
more complete picture of the nature of the 
dynamical interaction between the binary 
companion and the exoplanet. In addition 
to the system’s astrometry, their radial ve¬ 
locities should also be monitored. This 
may reveal a trend in the RV caused by 
the stellar companion. Knowledge of the 
RV trend coupled with measured astrome¬ 
try will enable the computation of the full 
orbital solution. 
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